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We study super-Planckian near-field heat exchanges for multilayer hyperbolic 
metamaterials using exact S-matrix calculations. We investigate heat exchanges 
between two multilayer hyperbolic metamaterial structures. We show that the super- 
Planckian emission of such metamaterials can either come from the presence of sur- 
face phonon-polaritons modes or from a continuum of hyperbolic modes depending 
on the choice of composite materials as well as the structural configuration. 

In the last few years several fascinating experiments have demonstrated that for small 
separation distances compared with the thermal wavelength the thermal radiation exchanged 
between two hot bodies out of thermal equilibrium increases dramatically compared with 
what we observe at large distances and can even exceed the well-known Stefan-Boltzmann 
law by orders of magnitude^-. Accordingly, thermal emission is in that case also called 
super-Planckian emission emphasizing the possibility to go beyond the classical black-body 
theory. There are several promising applications of super-Planckian emitters ranging from 
thermal imaging^- and thermal rectification/management)^— to near-field thermophoto- 
voltaics^~— . This has triggered many studies on the possibilities of tailoring and controlling 
the super-Planckian radiation spectrum by means of designing the material properties^— , 
using phase-change materials^ or 2D systems such as graphene^ 1 ^, for instance. 

Recently, it was shown that hyperbolic metamaterials can lead to broad-band photonic 
thermal conductance inside the material itself— and between two hyperbolic materials only 
separated by a vacuum gap2&. Further Nefedov et al. considered nanorod-like structures 
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made of nanotubes which are interlocked and highlighted a giant radiative heat flux which 
could be utilized for near-field thermophotovoltaics energy conversion^. Finally, Guo et 
al. have studied the energy density produced by the thermally fluctuating fields close to a 
hyperbolic structure and found a broadband near-field contribution from which they have 
concluded that super-Planckian emission will be broad-band for hyperbolic materials^. This 
is in accordance with the findings in Ref.— for the energy exchange between two hyperbolic 
nanowire structures. 

The aim of this letter is to show that the surface modes supported by the topmost layers of 
phonon-polaritonic metamaterials can give the dominant contribution to the super-Planckian 
emission. As was shown in Ref.— materials which have a broad hyperbolic frequency band 
as predicted from effective medium theory can support surface modes inside these frequency 
bands as well which will compete with the hyperbolic modes^i^. In particular, we will 
show that for the realization of a hyperbolic metamaterial as studied in Ref.— the main 
contribution to super-Planckian radiation is not necessarily due to hyperbolic modes but 
can be due to surface modes depending on the choice of the topmost layer. We will show 
that in order to allow for broad-band super-Planckian emission by hyperbolic modes, mainly, 
it is important to use a material for that topmost layer which does not support surface modes 
in the thermal freqeuency range. 

Before we start to study the super-Planckian thermal radiation let us first recall the 
concept of indefinite or hyperbolic materials^— . Such materials are first of all a special 
class of uni-axial anisotropic materials. For uni-axial materials the permittivity e\\ parallel 
to the optical axis is different from the permittivity e± perpendicular to the optical axis. 
For hyperbolic materials one can find frequency bands where en and e± have different signs, 
i.e. || < 0. Thus the dispersion relation of the photons in such a material^ 

K 2 k 2 I,! 2 

- + ^ = - (1) 
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describes a hyperbolic function rather than an ellipse as for usual anisotropic materials; 
here k (k z ) is the wavevector inside the hyperbolic medium perpendicular (parallel) to the 
optical axis which is assumed to point in ^-direction. Such metamaterials can for example 
be designed by multilayer structures, since in the long wavelength regime such structures 
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can be described as homogeneous anisotropic media with the effective permittivities 

e± = e 1 / + e 3 (l-/) J (2) 

where e% and e 2 are the permittivities of the two layer-materials and / is the filling fraction 
of the topmost material 1, i.e. / = h/(h + h)- The effective permittivities allow for a 
calculation of the hyperbolic frequency bands of the multilayer structure where e±e\\ < 0. 
There are in general two different kinds of bands: a frequency band Ai where e\\ < and 
e± > and a frequency band A 2 where eii > and e± < 0. As will become clear in the 
following the such calculated frequency bands Ai and A2 can also support surface modes, 
which are not taken into account in the effective description^. 




Figure 1: Sketch of the geometry of two hyperbolic multilayer materials separated by a vacuum 
gap. 

In order to study super-Planckian radiation we consider the geometry depicted in Fig. [TJ 
The heat transfer coefficient h(d) between the two metamaterials which are assumed to be 
at local thermal equilibrium can be determined by^ 

— f(co,T)J2 J-^T 3 ( Ut K ] d) = j( —f(u,T)H(L,,d) (4) 

where f(u,T) = (Hoo) 2 /(k B T 2 )e nw / ksT /(e hw / kBT - l) 2 . T s (uj,K;d) and %(uj,K;d) are the 
energy transmission coefficients for the s- and p-polarized modes which can be easily de- 
termined for semi-infinite materials, anisotropic materials and multilayer structures'^—. 
Here we use the standard S-matrix approach as in Refs.— & to calculate the amplitude re- 
flection coefficients Tj of our multilayer structures from which we can easily determine the 
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energy transmission coefficients 

f(i-y 2 )7IAf, K<ujjc 

Tj(u,K\d) = < (5) 
I 4[Im(r i )] 2 e- 2|fczol 7|D J f , k > u/c 

including the contributions of the propagating modes with k < u/c and the evanescent 
modes with k > u/c. Here Dj — 1 — rjrje 2lkz ° d is a Fabry-Perot-like denominator with 
k^Q = fcg — ko = u/c. 

Now, let us consider a concrete example of a hyperbolic structure which is composed by 
layers of polar materials. Because these structures can support surface phonon-polaritons as 
well, they are also called phonon-polaritonic hyperbolic structures. We choose to consider 
the structure in Ref.— which is made of layers of SiC and SiCV In general amorphous SiC>2 
supports surface modes in the infrared as well as SiC, but to get results which are comparable 
with the calculations done in Ref.— we assume that es;o 2 = 3.9 adding a vanishingly small 
absorption. The optical properties of SiC are taken from Ref.—. The layer thicknesses are 
(a) l\ = 50 nm for the SiC layers and I2 = 150 nm for the silica layers so that the filling 
fraction is f — 0.25 and (b) h = I2 = 100 nm so that / = 0.5. For our exact S-matrix 
calculations we use N = 50 where the last layer is a semi-infinite layer with the material 
properties of the topmost layer. The hyperbolic frequency bands calculated from Eqs. (J2J) 
and (J3J are (a) A 1 = 1.495 - 1.623 • 10 14 rad/s and A 2 = 1.778 - 1.826 • 10 14 rad/s and (b) 
Ai = 1.495 - 1.712 • 10 14 rad/s and A 2 = 1.712 - 1.827 ■ 10 14 rad/s. 

In order to see the structure of contributing modes we have plotted the transmission 
coefficient T p (u,K;d) in Fig. [2j The horizontal dashed white lines mark the hyperbolic 
bands as determined from effective medium theory^. The solid white lines are the borders 
of the Bloch bands as determined from Bloch mode dispersion relation for p polarization^ 

2\e 1 k z2 e 2 k zl J ^ 
+ cos(A; zl / 1 ) cos(k Z 2h), 

with the permittivities 6j (i = 1,2) of the two layer materials and the wavevector along 
the optical axis in z direction k z \ = \/k\ei — k 2 . Note that k z>B is the Bloch wavevector 
inside the multilayer structure which can be approximated by its homogenized version k z in 
Eq. (|T]) together with Eqs. d2J) and in the regime where the effective description is valid. 
Only inside these Bloch bands one can find modes which are propagating modes inside the 
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Figure 2: Transmission coefficient Tp(u),K;d) from Eq. (0) for both SiC-SiC>2 multilayer structures 
(a) l\ = 50 nm and \<i = 150 nm, and (b) l\ = h = 100 nm for the interplate distance d = 100 nm. 

hyperbolic material. It can be seen that there are also very dominant modes outside the 
Bloch bands contributing significantly to the energy transmission. These modes are the 
coupled surface modes of the topmost SiC layers of each hyperbolic material which means 
that they are evanescent modes inside and outside the hyperbolic structure. 

The respective contribution of the modes inside and outside the Bloch bands to the 
spectral heat transfer coefficient H(oj,d) is plotted in Fig. |3]for a distance of d = 100 nm. 
From that figure it becomes apparent that within the hyperbolic frequency bands one has 
quite large contributions steming from modes outside the Bloch bands which are mainly the 
coupled surface modes of the topmost layers. Hence, for the chosen structure the broadband 
super-Planckian radiation from the hyperbolic frequency band is not due to hyperbolic 
modes only. The relative contribution of surface modes and all the other modes is plotted 
in Fig. H] where we show the heat transfer coefficient as a function of distance. From that 
figure it becomes obvious that for distances about 100 nm and smaller the heat flux is 
dominated solely by the coupled surface modes of the topmost layers showing a typical 1/d 2 
dependences^. Whereas for larger distances the heat flux is dominated by the contributions 
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inside the Bloch bands. These contributions are on the one hand hyperbolic modes steming 
from frequencies inside the hyperbolic bands Ai and A 2 . On the other, for frequencies 
outside the frequency bands Ai and A 2 the modes are usual propagating or frustrated total 
internal reflection modes. Note that for distances of the order of max(/i, / 2 )/7r the Bloch- 
mode contribution reaches a maximum. This can be attributed to the large wavevector 
cutoff by the edge of the Bloch bands which can be understood as the inset of nonlocal 
effects since for such distances the main wavevector contributions to the thermal emission 
are of the order ir/d. 

To quantify the heat flux mediated by the hyperbolic modes we plot in Fig. [5] the different 
contributions of the modes inside and outside the Bloch bands, and the contribution from 
the hyperbolic modes separately. The separate contributions hs^d), tiNB(d), and hh m (d) to 
the heat transfer coefficient are normalized to the total heat transfer coefficient h tot (d) = 
hB(d) + hNB{d). Apparently, in both configurations the contribution of the hyperbolic modes 
is for all chosen distances smaller than 35%. This is a rather small value for a hyperbolic 
structure which is constructed for the purpose of enhancing the thermal radiation by the 
hyperbolic-mode contribution. 



1500 



g 1000 

500 



■total 
Bloch 

■Non-Bloch 



1.5 



1.6 1.7 

u (10 14 rad/s) 






2000 




1500 
















1000 






^3 






500 




9 



total 






(b) 


Bloch 








' Non-Bloch 










! A 2 

' 








— Tl 


\ 

. ^ v_ 





1.5 



1.6 1.7 

u (10 14 rad/s) 



1.9 



Figure 3: Spectral heat transfer coefficient H(u>, d) defined in Eq. @ normalized to the black-body 
result -Hbb(w) = u 2 /(2-kc 2 ) for both SiC-Si0 2 multilayer structures (a) l\ = 50 nm and Z 2 = 150 nm, 
and (b) h = h = 100 nm for the interplate distance d = 100 nm. Here we choose T = 300 K. The 
vertical dashed lines mark the borders of the hyperbolic frequency bands Ai and A 2 . 



Now, let us see if the dominant surface-mode contribution vanishes when choosing the 
passive Si0 2 -layer as the topmost layer. Here, it is important to keep in mind that Si0 2 
supports surface modes in the infrared. We assume here that it can be described by a 



7 




Figure 4: Heat transfer coefficient h{d) from Eq. Q as a function of interplatedistance d using 
T = 300 K for both SiC-Si02 multilayer structures (a) l\ = 50 nm and I2 = 150 nm, and (b) 
h = h = 100 nm. The heat transfer coefficient is normalized to the black-body value /ibb = 
6.1 Wm _2 K _1 . The contributions from the Bloch bands and from regions outside the Bloch bands 
are shown separately. 

constant permittivity in the frequency band of interest (it is in this sense "passive") in 
order to compare our results to existing results in the literature. Hence, we repeat the same 
calculations for the same structure as before but with the difference that for both hyperbolic 
structures the topmost layer is SiC>2 followed by SiC, etc. The results for the spectral heat 
transfer coefficient are shown in Fig. [6] (a). There is still a surface-mode contribution, but it is 
very small compared to the Bloch-mode contributions. Finally, from the distance dependent 
results in Fig. (b) and (c) it can be seen that the super-Planckian radiation is mainly due 
to Bloch modes, i.e. frustrated total internal reflection modes and hyperbolic modes. In 
particular, the contribution of the hyperbolic modes can be larger than 50% in the strong 
near-field regime for distances of about 10 nm. 

As is obvious from Fig. [6] (b) the overall heat flux is only one order of magnitude larger 
than that of a black body so that the hyperbolic material considered here and in Ref.— is 
a poor near-field emitter compared to the previous structures with SiC as topmost layer. 
But there is a simple method for increasing the hyperbolic contribution by just making the 
thickness of the layers smaller. Then, the border of the Bloch bands will shift to larger 
wavevectors which results in a broadband contribution to the transmission coefficient for 
larger wavevectors and hence to a larger thermal radiation. In Fig. [7] we show the heat 
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Figure 5: Heat transfer coefficients /ijvb> an d hh m of the Bloch modes, the modes outside the 
Bloch bands and the hyperbolic modes normalized to the total heat transfer coefficient htot(d) = 
hB(d)+fiNB(d) as a function of distance. Again we show both cases (a) l\ = 50 nm and I2 = f50nm, 
and (b) h = l 2 = 100 nm, and set T = 300 K. 

flux for hyperbolic structures with a filling factor of 0.5 but layer thicknesses of 100 nm, 
50 nm and 5 nm. It can be seen that the heat flux increases by orders of magnitude in 
the strong near-field regimes, i.e., for distances smaller than 100 nm, when making the 
layers thinner. We have checked that the main contribution is due to hyperbolic modes 
in that regime (not shown here). Further studies have to find an optimized design and 
optimal composite materials in order to further improve the thermal radiation properties 
of hyperbolic materials to attain thermal heat fluxes wich are as large as the heat flux by 
surface modes or even larger. Note that in Ref.— such a structure was proposed on the basis 
of an effective description. 

In conclusion, we have studied the super-Planckian emission of hyperbolic structures by 
using the fluctuational electrodynamics theory combined with the S-matrix method. It has 
been shown that to properly describe the energy exchanges it is of crucial importance not 
only to choose a good combination of material composites for having broad-band super- 
Planckian radiation but also to use a passive material as topmost layer, i.e. a material 
which does not support surface mode resonances within the thermally accessible spectrum. 
Also, we have shown for multilayer structures that the thickness of layers determines the 
wavevector cutoff of the Bloch band so that it appears clearly advantageous to use thin 
layers with elementary thicknesses I1J2 <§C d to observe a large super-Planckian emission 
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at a given distance d from the surface. These findings provide the basis for realizing an 
optimized design of hyperbolic thermal emitters with broad-band super-Planckian spectra. 
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Figure 6: (a) Spectral heat transfer coefficients H(lo, d) between two hyperbolic materials with 
Si02 as topmost layer choosing d = 100 nm. The vertical dashed lines mark the borders of the 
hyperbolic frequency bands Ai and A2. (b) Heat transfer coefficients h(d) for the same materials 
setting T = 300 K normalized to the black-body value /ibb — 6.1 Wm _2 K _1 . Finally in (c) we plot 
the relative contributions of the Bloch modes, Non-Bloch modes and the hyperbolic modes. 
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Figure 7: The heat transfer coefficient for the structure with the passive material as topmost layer 
for different layer thicknesses l\ = h (f = 0.5) of 100 nm, 50 nm, and 5nm normalized to the 
black-body value /ibb = 6.1 Wm _2 K _1 . 
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